First-principles investigation of graphene fluoride and graphane 
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Different stoichiometric configurations of graphane and graphene fluoride are investigated within 
density functional theory. Their structural and electronic properties are compared, and we indicate 
the similarities and differences among the various configurations. Large differences between graphane 
and graphene fluoride are found that are caused by the presence of charges on the fluorine atoms. A 
new configuration that is more stable than the boat configuration is predicted for graphene fluoride. 
We also perform GW calculations for the electronic band gap of both graphene derivatives. These 
band gaps and also the calculated Young's moduli are at variance with available experimental data. 
This might indicate that the experimental samples contain a large number of defects or are only 
partially covered with H or F. 
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I. INTRODUCTION 

Two-dimensional crystals have been given a large 
amount of attention since the isolati on o f one-atom-thick 
materials by Novosolov et al. in 2004P2] Graphene, a sin- 
gle layer of graphite, has attracted by far the most atten- 
tion because of the high crystal quality of the graphene 
samples and its fascinating electronic properties.^ These 
properties make it a promising candidate to use as a ba- 
sic material for future electronics applications!^ However, 
the use of graphene for applications in electronics suffers 
from a major drawback: graphene is, in its pristine state, 
a zero-bandgap semiconductor and this gapless state ap- 
pears to be rather robust. Several ways have been ex- 
plored to induce a finite band gap in graphene. It was 
found experimentally that a band gap can be opened by 
confining the electrons in nanoribbons^ or by a pplying a 
potential difference over a graphene bilayer!^ 

The chemical modification of gra phe ne is another 
promising way to create a band gapP*^ When radicals 
such as oxygen, hydrogen, or fluorine atoms are adsorbed 
on the graphene surface they form covalent bonds with 
the carbon atoms. These carbon atoms change their hy- 
bridization from sp 2 to sp 3 , which leads to the opening of 
a band gap (similar as in diamond). The adsorbed rad- 
icals can attach to the graphene layer in a r andom way, 
as is the case in graphene oxide (GO) JSE21 or they can 
form ordered patterns. In the last case, new graphene- 
based 2D crystals are formed with properties that can 
vary greatly from their parent material. This has been 
found to be the case for hydrogen and fluorine adsorbates. 
The new 2D crystals that are exp ected to form in those 
cases^ have been named graphane^^ and graphene flu- 
oride (or fluorographeneJPJJ, respectively. 

Following this rou te, m ulti-layer graphene fluoride was 
recently synthesized pESl anc j its structural and electronic 
properties were studied. A strongly insulating behav- 
ior was found with a room temperature resistance larger 
than lOGfi, which is consistent wit h the existence of a 
large band gap in this new material. 9 10 Only a partial 



fluorine coverage of the graphene multi-layer samples was 
achieved in these experiments. The F/C ratio was esti- 
mated to be 0.7 in Ref. [9] and 0.24 in Ref. 10, according 
to weight gain measurements. 

An important step forward in creating fully covered 
two-dimensional graphene fluoride samples was recently 
achieved in Ref. [TT] The obtained single-layer graphene 
fluoride exhibits a strong insulating behavior with a room 
temperature resistance larger than lTf2, a strong temper- 
ature stability up to 400 °C, and almost a complete dis- 
appearance of the graphene Raman peaks associated with 
regions that are not fully fluorinatedP^The graphene Ra- 
man peaks do not disappear completely, however, which 
could be an indication of the presence of defects in the 
sample, such as a small portion of carbon atoms not 
bonded to fluorine atoms. It was also found experimen- 
tally that fluorographene has a Young's modulus of rj 
100 N/m, and the optical measurements suggest a band 
gap of sa 3eV. 

In Ref. [H|it was demonstrated that single-side adsorp- 
tion is also possible and that it probably results in a crys- 
talline C4F structure with a large band gap. 

On the theoretical side, first-principles studies on 
graphene monofluoride started in 1993, motivated by 
available experiments on graphite monofluoride. Us- 
ing density functional theory (DFT) calculations, it was 
shown in Ref. [17]that the chair configuration of graphene 
fluoride is energetically more favorable than the boat con- 
figuration by 0.145 eV per CF unit (0.073 eV/atom), 
while a transition barrier of the order of 2.72 eV was 
found between both structures. Due to the small dif- 
ference in formation energy and the large energy barrier 
between both configurations, it was argued that the kine- 
matics of the intercalation could selectively determine the 
configuration, or that there could also be a mixing of both 
configurations in the available experiments. By using 
the local density approximation (LDA) for the exchange- 
correlation functional a direct band gap of 3.5 eV was cal- 
culated for the chair configuration in Ref. [17l However, 
it is well known that DFT generally underestimates the 
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band gap. Recent calculations used the more accurate 
GW approximation and found a much larger band gap of 
7.4 eV for the chair configuration of graphene monoflu- 
oride (Ref. ITS")) . This theoretical value is twice as large 
as the one obtained experimentally for graphene fluoride 
in Ref. [11] which is 3eV. The experimental value for 
the Young's modulus as found in Ref. [TT|(?s 100 N/m) is 
also half the value obtained recently from first-principles 
calculations in Ref. [T!|] (~ 228 N/m) for the chair config- 
uration of graphene fluoride. It is worth noting that the 
experimental and theoretical values of the Young's 
modulus of graphene only differ in a small percentage. 

Possible reasons for the disagreement between the ex- 
perimental values and the ab initio results for the Young's 
modulus and the band gap of graphene fluorine could be: 
i) the presence of a different configuration or a mixture of 
them in the experimental samples, or ii) the presence of 
defects, which could decrease the size of both the Young's 
modulus and the band gap from the expected theoretical 
values. 

In this paper, we investigate various possible crystal 
configurations for both graphene-based two-dimensional 
crystals, graphene fluoride and graphane, and we exam- 
ine their structural, electronic, and mechanical proper- 
ties. In the case of graphene fluorine, we found a new 
configuration not considered before that has a lower en- 
ergy than the boat configuration. This new configura- 
tion, which we call the zigzag configuration, is energeti- 
cally less favorable than the chair configuration by only 
0.073 eV per CF unit (0.036 eV/atom). We calculated 
the Young's modulus and the band gap (both with GGA 
and in the GW approximation) for the different config- 
urations. The disagreements between experimental and 
ab initio calculations for graphene fluoride persist inde- 
pendently of the considered configuration. These results 
imply that the available experimental samples probably 
contain a large number of defects, such as a portion of 
carbon atoms not bonded to fluorine atoms, that decrease 
the value of both the Young's modulus and the band gap 
from the expected theoretical values. 

The paper is organized as follows: first we describe the 
computational details of our first-principles calculations. 
Then we investigate the stability and structural prop- 
erties of the different configurations of both graphene 
derivatives. To conclude, the elastic and electronic prop- 
erties of the different structures are discussed. 



II. COMPUTATIONAL DETAILS 

We examine different graphane and graphene fluoride 
configurations with the use of ab initio calculations per- 
formed within the density functional theory (DFT) for- 
malism. The generalized gradient approximation (GGA) 
of Perdew, Burke, and Ernzcrhof (PBEJ^is used for the 
exchange-correlation functional and a plane wave basis 
set with a cutoff energy of 40 Hartree is applied. The 
sampling of the Brillouin zone is done with the equiva- 



lent of a 24 x 24 x 1 Monkhorst-Pack k-point gricP^l for 
a graphene unit cell and we use pseudopotentials of the 
Troullier-Martins typeP^ Since periodic boundary condi- 
tions are applied in all 3 dimensions a vacuum layer of 
20 Bohr is included to minimize the (artificial) interac- 
tion between adjacent layers. All the calculations were 
performed with the ABINIT codePS 

The reported quasiparticle corrections for the band 
gap are obtained using the YAMBO codeP^Here the first- 
order quasiparticle corrections are obtained using Hedin's 
GW approximation^ for the electron self-energy. Be- 
cause we are treating two-dimensional systems, the spu- 
rious Coulomb interaction between a layer and its im- 
ages should be avoided, as this causes serious convergence 
problems. Therefore we use a truncation of this interac- 
tion in a box layout, following the method of Rozzi et 
al.W* The remaining singularity is treated using a ran- 
dom integration method in the region near the gamma 
pointPSl Nevertheless, a larger separation between the 
layers is necessary, so a value of 60 Bohr is used for these 
calculations. 



III. RESULTS 

We studied four different stoichiometric configurations 
for both graphane and graphene fluoride in which every 
carbon atom is covalently bonded to an adsorbate in an 
equivalent way, i.e. every carbon/adsorbate pair has the 
same environment. These configurations are schemat- 
ically depicted in Fig. [I] and we will refer to them as 
the 'chair', 'boat', 'zigzag', and 'armchair' configuration. 
The chair and boat configurations have been well investi- 
gated before, but the zigzag and armchair configurations 
are rarely examined for graphane^ and we are not aware 
of any studies for fluorographene. The names of these last 
two configurations have been chosen for obvious reasons 
(see Fig. [ijc) and[ljd)). After relaxation, the different 
configurations appear greatly distorted when compared 
with the schematic pictures of Fig. [Tj so these figures 
should only be regarded as topologically correct (see Fig. 

A. Stability analysis 

To examine the stability of the different configurations, 
we make use of the formation energy of the structures 
and the binding energy between the graphene layer and 
the adsorbates. We define the formation energy, Ef, as 
the energy per atom of the hydrogenated or fluorinatcd 
graphene with respect to intrinsic graphene and the cor- 
responding diatomic molecules H 2 and F 2 . The binding 
energy, Eb, is defined with respect to graphene and the 
atomic energies of the adsorbates and is calculated per 
CH or CF pair. The results are summarized in Table [I] 

As has been reported before, the chair configuration is 
the most stable one for both graphane^ and graphene 
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FIG. 1: (Color online) Four different configurations of 
hydrogen/fluorine-graphene: a) chair, b) boat, c) zigzag, and 
d) armchair configuration. The different colors (shades) rep- 
resent adsorbates (H or F) above and below the graphene 
plane. The supercell used to calculate the elastic constants is 
indicated by the dashed box. 



TABLE I: The formation energy E/, the binding energy Ef,, 
and the relative binding energy AE/ (with respect to the 
most stable configuration) for different hydrogenated and flu- 
orinated graphene configurations. The energies are given in 

eV. 

chair boat zigzag armchair 

graphane 

E b -2.481 -2.378 -2.428 -2.353 

E f -0.097 -0.046 -0.071 -0.033 

AEf 0.000 0.051 0.027 0.064 

fluorographene 

E b -2.864 -2.715 -2.791 -2.673 

E f -0.808 -0.733 -0.772 -0.712 

AEf 0.000 0.075 0.036 0.095 



experiment allyP^ 



B. Structural properties 

Besides the large difference in formation energy there 
are also pronounced structural differences between both 
graphene derivatives. The structural parameters for the 
different configurations of graphane and fluorographene 
are shown in Table ITU Note that all the structures are de- 
scribed in an orthogonal supercell, as illustrated in Fig. 
[TJ for ease of comparison. The results for the chair config- 
uration agre e well w ith previous theoretical calculations 
for graphan e! 16 ! 30 * 31 ! and graphene fluorideP^ 

TABLE II: Structure parameters for the different hydro- 
genated and fiuorinated graphene derivatives. Distances are 
given in A and angles in degrees. The distance between neigh- 
boring C atoms, dec, an d the angles, 9ccx, ar e averaged over 
the supercell. 





chair 


boat 


zigzag 


armchair 


graphane 










a* /a/3 


2.539 


2.480 


2.203 


2.483 


a v /n y 


2.539 


2.520 


2.540 


2.270 


dcH 


1.104 


1.099 


1.099 


1.096 


dec 


1.536 


1.543 


1.539 


1.546 


#CCH 


107.4 


107.0 


106.8 


106.7 


&CCC 


111.5 


111.8 


112.0 


112.1 


fluorographene 










a x /V3 


2.600 


2.657 


2.415 


2.662 


Cly/fly 


2.600 


2.574 


2.625 


2.443 


dcF 


1.371 


1.365 


1.371 


1.365 


dec 


1.579 


1.600 


1.585 


1.605 


#CCF 


108.1 


106.0 


104.6 


104.2 


#ccc 


110.8 


112.8 


113.9 


114.2 



fluoride^. The zigzag configuration is found to be more 
stable than the boat and armchair configurations and its 
formation energy is only slightly higher than that of the 
chair configuration: for both graphene derivatives the 
difference in formation energy, AEf, between chair and 
zigzag is of the order of the thermal enenrgy at room tem- 
perature (26 meV). The energy differences between the 
various configurations are more pronounced for graphene 
fluorine than for graphane but they are of the same order 
of magnitude. 

When we compare graphane and fluorographene, the 
binding energy of hydrogen and fluorine appears to be 
rather similar (2.5 eV compared to 2.9 eV) but there 
is a huge difference in the formation energy (0.1 eV 
compared to 0.9 eV). This is a consequence of the large 
difference in the dissociation energy between hydrogen 
and fluorine molecules. The formation energy as defined 
above can be regarded as a measure of the stability 
against molecular desorption from the graphene surface. 
Therefore graphene fluoride is expected to be much 
more stable than graphane as has indeed been observed 



It is also useful to compare the interatomic distances 
and bond angles with those of graphene and diamond. 
Therefore we calculated these using the same formal- 
ism as described above (section [n| . The C-C bond has 
a length of 1.42 A for graphene compared to 1.54 A 
for diamond, and the bond angles are 120° and 109.5°, 
respectively. Notice that both graphane and fluoro- 
graphene resemble much closer the diamond structure 
than graphene, which is not surprising since the hy- 
bridization of the carbon atoms in these structures is the 
same as in diamond, i.e. sp 3 . The C-C bond length for 
the graphane configurations is similar to the one in dia- 
mond, but dec m fluorographene is about 0.05 A larger. 
This can be explained from a chemical point of view as 
due to a depopulation of the bonding orbitals between 
the carbon atoms. The depopulation of these bonding 
orbitals results from an electron transfer from the carbon 
to the fluorine atoms due to the difference in electroneg- 
ativity between C and F. We used a Hirshfeld-based 
methocPsHm to calculate this charge transfer and found 
it to be AQ ss 0.3e. The charge transfer in graphane is 
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much smaller because of the similarity between the elec- 
tron affinity of C and H. 
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FIG. 2: (Color online) Fluorographene in the zigzag (a) and 
armchair (c) configuration. The nearest neighbor bonds of one 
F atom are indicated with dotted lines to show the symmetry 
of the superlattices (as shown in (b) and (d)). b) The hexag- 
onal superlattice which is formed in case of chair and zigzag 
configurations, d) The cubic superlattice which is formed in 
case of boat and armchair configurations. 



The fact that the fluorine atoms are negatively charged 
has an appreciable influence on the structure of graphene 
fluoride when compared to graphane. This can, e.g., be 
seen from the sizes of the different bond angles. The 
bond angles (and also the bond lengths) in the chair con- 
figuration can be regarded as the ideal angles (lengths) 
for these structures because they can fully relax. The 
other configurations will try to adopt these ideal bond 
angles and it can be seen from Table [TT] that this is in- 
deed the case for the graphane configurations. The flu- 
orographene configurations, on the other hand, appear 
to be somewhat distorted because their bond angles are 
(relatively) far from ideal. This is probably caused by 
the repulsion between the different fluorine atoms as can 
be demonstrated when focusing only on the positions of 
the F atoms. The fluorine atoms appear to form hexago- 
nal or cubic superlattices depending on the configuration 
(see Figs.J^b) and^d)). This is trivial in the case of 
the chair (and maybe the boat) configuration but not so 
for the others (see Figs. [2ja) and[2|c)). These superlat- 
tices are not perfect (deviations of a few percent), but 
are much more pronounced than in the case of graphane. 

So it seems that, at the cost of deforming the bond- 
ing angles, F superlattices are formed to minimize the 
electrostatic repulsion between the charged F atoms. 



C. Elastic strain 

Graphene and its derivatives graphane and fluoro- 
graphene can be isolated and made into free-hanging 
membranes. This makes it possible to measure the elas- 
tic constants of these materials from nanoind entation ex- 
periments using an atomic force microscope! 8 ! 11 ^ The 
experimental elastic constants can be compared to first- 
principles calculations which gives us information about 
the purity and structural crystallinity of the experimen- 
tal samples. Therefore we calculated the (2D) Young's 
modulus, E' , and the Poisson's ratio, v, of the differ- 
ent graphane and fluorographene configurations, which 



we list in Table III The Young's modulus and the Pois- 
son's ratio of graphene are found to be E' = 336 Nm" 1 
and v = 0.17, respectively, which corresponds well to the 
experimental value, E' exp — 340 ± 50 Nm -1 , and other 
theoretical results JEJI25I 



TABLE III: Elastic constants of the different hydrogenated 
and fluorinated graphene derivatives. The 2D Young's mod- 
ulus, E' , and Poisson's ratio, v, are given along the cartesian 
axes. E' is expressed in Nm - . 





chair 


boat 


zigzag 


armchair 


graphane 










E'c 


243 


230 


117 


247 


K 


243 


262 


271 


142 


v x 


0.07 


-0.01 


0.05 


-0.05 


Vy 


0.07 


-0.01 


0.11 


-0.03 


fluorographene 










E' c 


226 


238 


240 


215 


K 


226 


240 


222 


253 




0.10 


0.00 


0.09 


0.02 


Vy 


0.10 


0.00 


0.11 


0.02 



The 2D Young's modulus of graphane and fluoro- 
graphene is smaller than that of graphene. The E' of the 
chair and boat configurations of both graphene deriva- 
tives are about 2/3 the value of graphene which makes 
them very strong materials. The Young's modulus for 
the zigzag and armchair configurations of graphane are 
highly anisotropic with values that are roughly halved 
along the direction that shows the largest crumpling (see 
Figs. [2] (a) and[2jc)). The situation is completely dif- 
ferent for fluorographene where the Young's modulus is 
more isotropic. The values that are found for the chair 
configurations agree well with recent calculations (245 
Nm -1 and 228 Nm -1 for graphane and fluorographene, 
respectively, in Ref. I19p. Nair et al. performed a nanoin- 
dentation experiment on fluorographene^ and measured 
a value of 100 ± 30 Nm -1 for E' FG . This value is approxi- 
mately half the theoretical value, which suggests that the 
experimental samples contain a lot of defectsP^ 

The Poisson's ratio shows a similar behavior as the 
Young's modulus. The knowledge of E' and v allows 
us to calculate all the other 2D elastic constants^ such 
as the bulk, K' — E' /2(l — v), and shear modulus, 
G' = E' /2(l + v). For the chair configurations we find 
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Kft G = 131 Nm" 1 
and = 126 

graphene fluoride. 



and G[ 



HG 



114 Nm" 



Nm" 1 and G 1 , 



FG 



1 for graphane, 
103 Nm" 1 for 



D. Electronic properties 

Graphene is a zero-gap semiconductor but its deriva- 
tives, like graphane and fluorographene, have large band 
gaps, similar to diamond. In Table |IV[ the band gaps 
of the configurations under study are given. We also 
performed GW calculations because GGA is known to 
underestimate the band gap. The GGA and GW results 
show different behavior for the variation of the band gap 
among the different configurations. Note that this indi- 
cates that it is not straightforward to deduce qualitative 
trends from GGA as is often done in the literature. But, 
overall, we may conclude that the band gap is more or 
less independent of the configuration and that its size is 
roughly twice as large for GW as compared to GGA. 

The GGA results give a band gap of 3.2 eV for the 
most stable fluorographene configuration which is in ac- 
cordance with the experimental result of ~3 eV as found 
in Ref. 11. However, this value is much smaller than the 
(more accurate) GW band gap of 7.4 eV, so that the the- 
oretical and experimental results differ by about a factor 
of two. This conflict might be resolved if the experimen- 
tal value is ascribed to midgap states due to defects in 
the system, such as missing H/F atoms (similar to what 
has been predicted for defected graphane 3 ^). 



TABLE IV: Electronic properties of the different configura- 
tions of hydrogenated and fluorinated graphene derivatives. 
The electronic band gap, E gap , is given for GGA and GW 
calculations. The ionization energy (IE) is also calculated. 
All the energies are given in eV. 





chair 


boat 


zigzag 


armchair 


graphane 










E gap (GGA) 


3.70 


3.61 


3.58 


3.61 


Egap (GW) 


6.05 


5.71 


5.75 


5.78 


IE (GGA) 


4.73 


4.58 


5.30 


4.65 


fluorographene 










E gap (GGA) 


3.20 


3.23 


3.59 


4.23 


Egap (GW) 


7.42 


7.32 


7.28 


7.98 


IE (GGA) 


7.69 


7.64 


7.85 


8.27 



The electronic band structure and the corresponding 
density of states of graphane and fluorographene in the 
chair configuration are shown in Fig. [3] Both band struc- 
tures look similar but there are also some clear differ- 
ences. In the case of fluorographene the parabolic band at 
the T-point, corresponding to quasi-free electron states, 
is at much higher energies which indicates a larger ion- 
ization energy for fluorinated graphene. This ionization 
energy (IE) is defined as the difference between the vac- 
uum level and the valence band maximum and an explicit 
calculation of this energy indicates a difference of about 
3 eV between graphane and fluorographene (see Table 



IV). This is a consequence of the negative charges on 
the fluorine atoms in fluorographene. We can also com- 
pare the IE values with the work function of graphene 
which is the same as its ionization potential (because 
graphene has no band gap) and has a value of 4.22 eV 
from GGA (this is somewhat smaller that the experimen- 
tal value 38 of 4.57±0.05 eV). It can be seen from Table [fv 
that the ionization energies of both graphene derivatives 
are higher than that of graphene (« 0.5 eV and 3.5 eV, 
respectively) , although the ionization energy of graphane 
is rather similar to graphene. 




r DOS 



FIG. 3: (Color online) The electronic band structure and the 
corresponding density of states (GGA) for the chair configu- 
ration of graphane (a) and fluorographene (b). The valence 
band maximum has been used as the origin of the energy 
scale. 



IV. SUMMARY AND CONCLUSIONS 

We investigated different configurations of the 
graphene derivatives fluorographene and graphane. The 
chair configuration is the most stable one in both cases, 
but the zigzag configuration has only a slightly higher 
formation energy and is more stable than the much more 
studied boat configuration. Fluorographene is found to 
be much more stable than graphane which is mainly due 
to a much higher desorption energy for F2 as compared 
to H2. We also demonstrated that there are structural 
and electronic differences that are caused by the charged 
state of the F atoms in fluorographene. 

When our results are compared to available experimen- 
tal data for fluorographene some discrepancies can be 
noticed: for all the configurations studied we find much 
larger band gaps in the electronic band structure and the 
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calculated Young's modules is much larger. This might 
indicate that the experimental samples still contain ap- 
preciable amounts of defects. The nature of these defects 
requires further investigation, but one can speculate that 
these defects consist of missing adsorbates, partial H/F 
coverage, or mixed configurations. 
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